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Isospin violation in the vector form factors of the nucleon 
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A quantitative understanding of isospin violation is an increasingly important ingredient for the 
extraction of the nucleon's strange vector form factors from experimental data. We calculate the 
isospin violating electric and magnetic form factors in chiral perturbation theory to leading and next- 
to-leading order respectively, and we extract the low-energy constants from resonance saturation. 
Uncertainties are dominated largely by limitations in the current knowledge of some vector meson 
couplings. The resulting bounds on isospin violation are sufficiently precise to be of value to on-going 
experimental studies of the strange form factors. 

PACS numbers: 12.39.Fe, 13.40.Gp, 14.20.Dh 



I. INTRODUCTION 

Since a nucleon has no valence strange quarks, the 
strange form factors provide information about the dy- 
namics of virtual quarks within a nucleon. In recent years 
the strange electric and magnetic form factors have been 
studied experimentally by the SAMPLE 0, A4 HQ, 
HAPPEX [113, and GO @ Collaborations^nd theo- 
retically via lattice QCD simulations B. la l p EH EH D3 > 
chiral perturbation theory lla. X\M JTa, Ila . fl7| . and hadron 
models [H OS HI SJ ^,^M^^(Fm a discussion 
of earlier theoretical studies, see Ref. j2||.) A central 
conclusion of recent research is that the strange electric 
and magnetic form factors are small, perhaps even small 
enough to require a careful understanding of the compet- 
ing effects from isospin violation. 

The inequality of up and down quarks, in terms of both 
mass and electric charge, produces effects that mimic the 
strange form factors. Extraction of the authentic strange 
quark effects from experimental data requires this isospin 
violation to be understood theoretically. There have been 
a number of theoretical studies of isospin violation in this 
context |28l |2l| y(J, |3j , but it is not yet clear that 
uncertainties in these results are negligible with respect 
to the strange form factors. Refs. |27l |28j used quark 
model discussions to determine the isospin violating con- 
tributions, and found them to be less than 1%; in fact the 
strange magnetic form factor received no contribution at 
all for vanishing momentum transfer. This null contribu- 
tion at Q 2 = is not due to any symmetry of QCD, and 
therefore the chiral perturbation theory (ChPT) study in 
Ref. [2I]] did find a non-zero contribution, but it includes 
an undetermined low-energy constant. One could think 
of extracting the physics of this constant, for example, 
from the light-cone baryon model of Ref. [3(| or from the 
quark model approach of Ref. [32, though it is difficult 
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to see how to match these models to the consistent low- 
energy effective theory description of ChPT in order to 
avoid double counting problems. 

In the present work we address two main goals. First 
we revisit the calculation of chiral loops, this time using 
both a manifestly covariant formalism and heavy- 
baryon ChPT. We find that the two formalisms offer 
complementary advantages for various aspects of this 
calculation. Since the isospin violating effects begin at 
rather high orders in the chiral expansion, leading order 
(LO) = 0(p 4 ) and next-to-leading order (NLO) = 0(p 5 ) 
(where p collectively stands for small parameters like mo- 
menta, the pion mass, or the electric coupling), the cross- 
check of the entire calculation in both formalisms is valu- 
able, and we are able to identify and correct some errors 
in Ref. 0. 

Our second major goal in the present work is to esti- 
mate the lone combination of low-energy constants that 
appears in the chiral perturbation theory expressions. 
For this, we employ resonance saturation. The physics 
producing isospin violation in the nucleon's vector form 
factors is seen to be p-u: mixing, and numerical values for 
the required vector meson couplings are obtained from 
recent dispersive analyses. 

A question that cannot be fully answered by our work 
relates to the convergence of the chiral expansion. With 
only LO and NLO results, it is impossible to make any 
definitive statement about convergence. A completely 
consistent extension beyond NLO seems unwarranted, 
since it would involve two-loop integrals and would in- 
troduce additional unknown low-energy constants. How- 
ever, given the apparent smallness of the strange form 
factors and of the isospin violating form factors, we can 
neglect contributions that are simultaneously isospin vi- 
olating and strange. This allows the discussion of isospin 
violation to occur within two-flavor ChPT, where con- 
vergence properties are known to be dramatically better 
than for the three-flavor theory. 

The present work thus provides a complete discussion 
of the physics at LO for the electric and at NLO for the 
magnetic isospin violating form factors, producing nu- 
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merical values for isospin violation in the magnetic mo- 
ment, electric radius, and magnetic radius, which are of 
direct relevance to the ongoing experimental studies of 
strange form factors. 

The outline of our presentation is as follows. We de- 
fine the form factors under consideration and their lead- 
ing moments in Section ^ The chiral representation of 
the Dirac form factors to leading order and of the Pauli 
form factors to next-to-leading order is (re-)derived in 
some detail in Section lTm In order to provide resonance- 
saturation estimates for unknown chiral low-energy con- 
stants, we calculate the contributions of p-uj mixing in 
Section lTvl which also allow for some insight into higher- 
order contributions beyond what is strictly calculated in 
ChPT. The pertinent formalism to do this consistently 
with chiral constraints is included. We conclude our find- 
ings in Section Some technical details are relegated 
to the appendices. 



II. BASIC DEFINITIONS 



Finally, the proton's neutral weak form factors G P E ^ M de- 
pend on one specific linear combination of isospin break- 
ing form factors, 



(l- ism 2 6 W )G E/M (t)-G E/M (t) 
~ G S E/M(t) — G E / M (t) , 



(5) 



where 9w is the weak mixing angle, G s E , M (t) are the 
strange vector form factors, and 



G 



E/AI\ 



(t) - G E/M (t) - G E/M {t) 



E/AI i 



(0) 



We will use an analogous notation also for the lead- 



ing moments, i.e. K u,d 



etc. Eq. (0) demon- 



strates that the isospin breaking form factors "simulate" 
strangeness form factors even in a world without strange 
quarks 0| , such that they have a direct impact on 
the accuracy to which the strangeness form factors can 
be extracted from data. 



We define the isospin symmetry breaking Dirac and 
Pauli form factors according to 

(p(p') + n(p') | -(u^u - dj^d) | p(p) + n(p)) 



= u(p>)\^F?(t) + l -^Ft(t)}u(p), 
L Ztun - 1 

(p(p') - n(p') | i (ti7 M u + dj^d) | p(p) + n(p)) 



(1) 



u(p') ^F{(t) 



u(p) , 



2uin 

}}, $ refer to isospin 



where g M = p^-p^, t = q 
breaking pieces in the isovector and isoscalar vector cur- 
rents, respectively. Our notation is linked to that used 
in Refs. by 

Ff{t) = ^- d Ff +n (t) , F({t) = u + d Ff- n {t) . 

The Sachs form factors are defined in terms of these as 

' -/-.;■"(/! . 

(2) 



G*/(t)=F^(t) + 



im\ 2 



N 



Current conservation implies if (0) = 0, while the Pauli 
form factors are normalized to the isospin breaking pieces 
of the (anomalous) magnetic moments, 



III. CHIRAL PERTURBATION THEORY 
A. Effective Lagrangians 

We have performed the relevant loop calculations both 
in the heavy-baryon formalism |35ll36| and in the infrared 
regularization scheme of relativistic baryon ChPT |33[ . 
For reasons of briefness, we only display the relativistic 
forms of the relevant Lagrangian terms in this section. 
The equivalent heavy-baryon forms can be found in the 
quoted references. 

The Goldstone boson (pion) Lagrangian is given, at 
leading order, by 



(u^ + X +) 



(7) 



= U collect the pion fields in 
iu'V^Uu^ with the covariant 



u'x u 



ux u contains 



7T7T 

where the matrices u 2 
the usual manner, -- 
derivative V M . The field x+ 
the quark mass terms, x = 2-B diag(TO u , m<j) + . . ., where 
B is linked to the quark condensate in the chiral limit, 
and F can be identified with the pion decay constant, 
F = Fn = 92.4 MeV. (. . .) denotes the trace in flavor 
space. 

A striking observation is the fact that the leading 
term generating the charged-to-neutral pion mass differ- 
ence l37l. 



if' (Q) 



(3) 



To complete the definition of the leading moments of 
these form factors, we define radius-like terms as the co- 
efficients of form factor terms linear in t, 



F: 



1/2 W 
E/Af ' 



F^(0) + pi'/ 2 t + O(t 2 ) 



G^ M (t) = G^/ M (0) + tf>f M t + 0{t 2 ) 



(4) 



4 2 i 7 , = ZF\Q 2 + 



M 2 + ~ M 2 a = 2Ze 2 F 2 



where Q± = ^(uQu^ ± u^Qu), Q — diag(e,0), does not 
feature in the following calculation: in nuclear physics 
language, this operator only breaks charge indepen- 
dence (independence under completely general rotations 
in isospin space), but not charge symmetry (rotations by 
7r about a fixed axis in isospin space, resulting in the 



3 



simultaneous exchange of u and d quarks, as well as pro- 
tons and neutrons) as required for the form factors in 
Eq. , so the pion mass difference alone will never gen- 
erate contributions to the form factors under investiga- 
tion, and we will neglect terms that are of second order 
in isospin breaking. 

The leading-order pion-nucleon Lagrangian is given by 

=#{^- TO + ^75}* , (8) 

where is the usual covariant derivative acting on the 
nucleon, m is the nucleon mass in the chiral limit, and 
gA can be identified with the axial coupling constant, 
gA = 1.26. Of the second-order pion-nucleon Lagrangian 
(that also contains all possible effects of virtual photons), 
we only quote the terms relevant for our calculation: 

C%. =*{c 5 x + + f 2 F 2 (Q + )Q + 

, (9) 

Here we have introduced the notation O = O — \{0) 
for traceless operators. F+, = iv F^u + uF^u^ contains 
the usual electromagnetic field strength tensor. There are 
two types of terms in Eq. (|5J): cq and c-j are related to 
the isovector and isoscalar anomalous magnetic moments 
(in the chiral limit) according to 

K V — K p — K n — Cq , K S = K p + K n = Cq + 2c-j , (10) 

with the experimental values n v = 3.706, k s = —0.120, 
while the operators proportional to C5 and f% generate 
the leading proton-neutron mass difference 

Am = m n -m p = 4B(m u - m d )c 5 + e 2 F 2 f 2 ■ (11) 

As all our results will be expressible in terms of Am, we 
do not have to care about the precise separation of this 
mass difference into its strong (oc C5) and electromagnetic 
(oc f 2 ) parts. We have neglected to spell out terms in 
Eq. © that lead only to a common nucleon mass shift 
and can be taken care of by replacing the nucleon mass 
in the chiral limit by its physical value, m — ► mjy = 
i(m p + m n ). 

Finally, the following fourth-order terms [SfiL l3^ | con- 
tribute to the isospin breaking form factors: 

4V = *{- f f I Wx + - f f(^x + ) 

+ h i0 F 2 (Q + )(Q + F+) (12) 
+ h^F 2 {Q + )Q + {F+ y )Y^. 

The terms proportional to &xot, contribute to «/, the 
terms multiplied by eios, ^40 to kP . Therefore, both 
form factors contain one counterterm each proportional 
to m u — md and e 2 . The constants eios, hio contain a 
divergent piece, 

6108 = e l08M + fle 10fs L , = /140(A) + f3h i0 L , 



where L has a pole in (d — 4), 

L = l6^{rf 1 4 + ^^- 1 - l0g4 ^}' 

A is the dimensional regularization scale, and je is Eu- 
ler's constant. The scale-dependent finite constants e\ os , 
h\ obey a renormalization group equation according to 

X dx e ™ w -~TM S ' X dx h4o{X) -~i6^- 

In order to render the isospin breaking magnetic mo- 
ments finite and scale- independent, we find the beta func- 
tions 

r = 2 9a c 5 o 9 2 A h 

In contrast, eio7 and /143 are finite and scale- independent. 

It should be mentioned that Ref. [23 included the 
A(1232) as an explicit field in ChPT. When the corre- 
sponding loop diagrams were computed, their momentum 
dependences were found to be negligible relative to the 
final error bands on the form factors. Since these effects 
would also be insignificant relative to the error bands 
that will be obtained in the present work, we omit the 
explicit A(1232) field from the outset. In the framework 
used here, its effects would show up via higher-order low- 
energy constants that are beyond the accuracy we aim at 
in this work. 



B. Power counting 

We repeat here some power counting arguments that 
were mostly already put forward in Ref. 29] . 

It is easy to see that the usual (isospin conserving) 
vector form factors receive polynomial contributions to 
the leading moments, i.e. magnetic moment, electric (or 
Dirac) radius, magnetic (or Pauli) radius, at chiral or- 
ders p 2 , p 3 , p 4 , respectively. Pion loop contributions 
start at 0(p 3 ), therefore such loop effects can generate a 
parameter-free leading-order prediction for the magnetic 
radius. The leading magnetic radius term has to scale 
like M" 1 ; this infrared singularity has been well-known 
for a long time 01 ■ 

From the previous section, it is obvious that polyno- 
mial contributions to the leading moments of the isospin 
violating form factors always have to include factors of ei- 
ther m u — nid or e 2 , therefore they are suppressed by two 
orders in the chiral power counting with respect to the 
isospin conserving ones. We therefore find precisely the 
terms of » in Eq. (|T2|l contributing to the magnetic 
moments, while operators of electric and magnetic radius 
type would arise at orders p 5 and p 6 . On the other hand, 
the pion loop diagrams with additional mass insertions of 
the C5, f 2 operators in Eq. 10, shown in Fig. start to 
contribute only one order higher than the same diagrams 




(a) (b) (c) 
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FIG. 1: One-pion-loop diagrams. Full/dashed/wiggly lines 
denote nucleons/pions/vector currents, respectively. The 
filled circles in diagrams (d) and (e) represent the magnetic 
couplings from * . The crossed diagrams of (c) and (e) 

are not depicted separately. Mass insertions on the nucleon 
propagator lines, yielding the physical proton and neutron 
masses, are not explicitly shown. 




FIG. 2: A pair of two-loop diagrams. The closed wiggly lines 
denote virtual photons. The sum of both contributions to the 
magnetic form factor drops out. 

without mass corrections, therefore leading loop contri- 
butions are also of order p 4 . This means that, while a 
prediction of the isospin violating magnetic moments is 
hampered by the a priori unknown „. counterterms, 
both electric and magnetic radii can be unambiguously 
predicted at leading order, and for the magnetic radii 
even the next-to-leading-order corrections are free of un- 
known parameters. The leading infrared singularity in 
the magnetic radius scales like Am/M^, the leading elec- 
tric and subleading magnetic radius terms like Am/M^. 

A calculation of the isospin violating Pauli (or mag- 
netic) form factor up to 0(p 5 ) is massively facilitated by 
the fact that no two-loop diagrams contribute. This can 
be seen as follows: 

1 . As mentioned in Section flll Al the pion mass differ- 
ence alone cannot generate charge symmetry break- 
ing terms. Therefore, diagrams with two pion loops 
would require another (subleading) charge symme- 
try breaking vertex or the nucleon mass difference 
in order to contribute, which would then, however, 
be at least of order p 6 . 

2. Two photon loops are of second order in isospin 
breaking and can be disregarded. In addition, it is 
easily seen in the heavy-baryon formalism that such 
diagrams with only leading-order photon couplings 
cannot generate the spin operators necessary for a 
magnetic contribution. 



FIG. 3: One-photon-loop diagrams with magnetic couplings. 
The sum of both diagrams vanishes. (The addition of a 
crossed right diagram is implied.) 



3. Finally, there might be diagrams with one pion and 
one photon loop. The only diagrams that generate 
a magnetic structure at 0(p 5 ) are of the type (a) in 
Fig. with one additional photon loop attached. 
However, it can be checked in the heavy-baryon 
formulation that the sum of the two diagrams in 
Figgis proportional to the anticommutator of the 
two Pauli-Lubanski spin operators stemming from 
the pion-nucleon couplings, and therefore it again 
only yields a contribution to the electric form fac- 
tor. The same mechanism can be checked for all 
other possible diagrams. 

Furthermore, also one-loop diagrams with isospin- 
breaking vertices (other than the nucleon mass difference 
insertion) can be ruled out of consideration: 

4. Tadpole graphs with isospin breaking couplings fail 
to generate infrared singularities proportional to 

M-\ 

5. The third-order pion-nucleon Lagrangian contains 
isospin breaking pion-nucleon coupling constants. 
However, these affect only the ir°NN coupling, 
hence they do not contribute in a type (a) diagram, 
while the remaining diagrams in Fig.QJare sublead- 
ing in their contributions to the Pauli form factor 
and therefore can only play a role at 0(p 6 ). 

Finally, the only one-photon loops contributing to the 
magnetic form factor are the ones depicted in Fig. [3] 
They have no i-dependence up to the chiral order con- 
sidered here, and it is easily calculated that their contri- 
bution to ^(0) exactly cancels. 

Therefore we have proven that the only infrared singu- 
lar contributions to the charge symmetry breaking form 
factor radii, and all contributions up to 0(p 4 ) for the 
Dirac and up to 0(p 5 ) for Pauli form factor in addition 

to the C^N-t* counterterms in Eq. are given by nu- 
cleon mass difference effects in the diagrams in Fig. ^ 



C. Chiral representation of the form factors 

In this section, we write down the chiral representa- 
tions of the charge-symmetry breaking form factors, the 
Dirac form factors to leading order, the Pauli form fac- 
tors up to next-to-leading order. We decompose all form 
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factors according to 

where current conservation dictates 



(13) 



Ff (0) = F( (0) = . 

For convenience, we define an overall dimcnsionlcss pre- 
factor 



C 



g\m^Am 
(4tt^)2 



1.4 x 10" 



(14) 



Then the ^-independent terms of the Pauli form factors 
are given by 



*?(0) 



■ AC 
+ 8mjv 
itf (0) = -AC 
+ 8m n 



1 + 21og 



AU tt(k v + 6)M W 
A 2m,N 

- 2B(m u - m d )e r 10S 

tt(k s + l)M n 



e 2 F 2 h r A0 



(15) 



m N 

e 2 F 2 h iz - 2B(m u - m d )e w7 



Altogether, we find for the ultimately required combina- 
tion 



G U M d (0) 



^(o)-f|(o) 

-AC 



tt(2k s - k v - A)M V 



2tojv 



CT > 



(16) 



where 



SmjV e 2 i ;l2 (/l43-/j4o) _2 - B ( m " _m 'i)( e 107-e^o8) 

(17) 

Note that in order to ease notation, we have suppressed 
the scale dependence of h\ Q , e\ 08) Kq t in Eqs. (|13|l - l(T7|l 
that is necessary to compensate for the chiral logarithms. 

Up to the order we are considering, the chiral rep- 
resentations of the form factors F^ 2 {t) display no t- 
dependence, 



F({t) = Fiit) = 



(18) 



This is due to the fact that only diagram (a) in Fig. 
contains a two-pion-cut and therefore generates momen- 
tum dependence in the low-energy region, and that the 
isoscalar vector current does not couple to pions at this 
order. The i-dependence of the form factors Ff(t), on 
the other hand, is given unambiguously in terms of pion 
loop contributions and can be written as 



Hit) = 



2g\M 7T Am 

^ 7T 

Ag 2 A mN Am 
F3 



70 W " 473 (*) 



M 

x[m~— 7o(*)-573(*) 
mN 



(19) 



(20) 



The explicit representations of the loop functions 
7o/ 3 (t) are given in Appendix IA II We note that the 
corresponding expressions reported in Ref. |29j | contain 
errors: Eqs. (53) and (54) in that work are too large 
by factors of 2 and 4 respectively; Eqs. 119fl and l|20() 
of the present work are the correct expressions, and we 
have obtained them separately from heavy-baryon ChPT 
and from infrared regularization. Expanding the loop 
functions to linear order in t, we find the radius term 
according to 



Pm 



2C_ 
3M* 
2C_ 



Am N 



m N 



Am 



N 



(21) 



where the first correction stems from F^(t) and the sec- 
ond from the subleading terms in F^it). Numerically, 
we find that, due to the large enhancement factor of 
77r/4 « 5.5, the next-to-leading-order correction reduces 
the leading prediction for the magnetic radius term by as 
much as 80%. 



IV. CONTRIBUTIONS FROM p-uj MIXING 

The missing ingredient in order to make the chiral 
representations of the previous section predictive is a 
value for the counterterm combination K^j T that is uncon- 
strained from symmetry arguments. As no phenomeno- 
logical information is available to fix these constants, we 
have to resort to model estimates in order to get a handle 
on them. The obvious model to use seems to be resonance 
saturation. 

The method to estimate chiral low-energy constants 
by including heavier resonance fields in the theory and 
matching their low-energy limit to higher-order opera- 
tors was first introduced systematically in the context 
of meson ChPT in Refs. [H EJ El It was shown to 
work rather well for the constants in c!£\ , in particular, 
those constants that receive contributions from vector 
and axial vector resonances were shown to be numerically 
dominated by these (a modern version of the concept of 
"vector meson dominance" ) . 

In Ref. El this method was extended to the pion- 
nucleon sector, nucleon resonances (in particular the 
A(1232)) were included, and it was demonstrated that 

(2) 

most of the constants in C]l n can also be well understood 
this way. However little is known about how well this 
procedure works for higher-order pion-nucleon coupling 
constants. In Ref. 01 it was found for the nucleon vec- 
tor form factors that in particular the isovector couplings 
are not well saturated by the p contribution alone. This 
observation agrees with the results of various dispersive 
analyses of these form factors (see e.g. Refs. |45l \A6L l47j 1 
that always require several resonances per channel in or- 
der to describe the data adequately. 
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Despite this rather pessimistic view, there is reason to 
believe that resonance saturation for the isospin violat- 
ing vector form factors might in principle work better 
than for the isospin conserving ones; we will comment on 
this point again at the end of Section IIVBI The only 
resonances that can contribute at tree level are vector 
mesons, and isospin violation is provided by the well- 
known mixing of the isovector p a and the isoscalar u> 
state. We lay out the formalism for incorporating vec- 
tor mesons, their mixing and coupling to vector currents 
and nucleons in the following section before presenting 
analytical and numerical results for the isospin violat- 
ing magnetic moments. Contributions to the higher mo- 
ments (electric and magnetic radii) will come as a benefit 
that allows us to estimate higher-order corrections to the 
leading chiral loop predictions. 



A. Formalism, Lagrangians 

We here discuss the necessary ingredients to describe 
the interaction of vector mesons with external currents as 
well as nucleons. We make use of the antisymmetric ten- 
sor formulation that was shown to be particularly useful 
in the context of ChPT and resonance saturation investi- 
gations [321 0| . Most of the literature is concerned with 
vector mesons in the context of SU(3) ChPT; we will here 
only be concerned with the SU(2) subsystem and adapt 
the formalism accordingly. In SU(2), the antisymmetric 
tensor V pv is given in terms of p and to fields according 
to 

The free Lagrangian for V^ v then takes the form 



(22) 



In particular, Eq. (|22(l leads to a common vector meson 
mass M p = M u = My. For numerical results, we will 
use the mass of the p, My = 770 MeV. The coupling of 
the vector mesons to external vector (and axial vector) 
sources is given by 



Cvj — 



Fy 
2V2 



(23) 



where T pv contains the electromagnetic field strength 
tensor proportional to the quark charge matrix Q = 
e/3 diag(2, — 1) (as opposed to the nucleon charge ma- 
trix Q used elsewhere in this text). Eq. correctly 
reproduces the SU(3) relation for the vector meson de- 
cay constants F p = 3F U = Fy- Phenomenologically, 
this relation is rather well fulfilled, F p = 152.5 MeV, 
F^ = 45.7 MeV. 

Mixing of p° and lo has been discussed in a formal- 
ism similar to the one presented here in Ref. |4jj|, albeit 



in the framework of SU(3). Here, we perform an anal- 
ogous construction for SU(2). We are not interested in 
(isospin-conserving) quark mass renormalization effects 
of the vector meson masses, but only in terms contribut- 
ing to p-uj mixing. One single such term can be con- 
structed, 



r — 



(24) 



where we have once more made use of the notation 
= V^v - ^(V^). If we match Eq. to the 

SU(3) result in Ref. [51] > a l so invoking quark counting 
rules |49j |. the mixing parameter can be identified with 
Vmix = 2My/B. In addition, there is a p-u transition 
through an intermediate photon via the interaction term 
Eq. 1)23(1 . Combining both effects, one finds the on-shell 
mixing amplitude 



Q puJ = 2M v (m u - m d ) + 



e 2 F v 



(25) 



The latest analysis of experimental data [5(j yields 

9 pw = (-3.75 ± 0.36) x 10~ 3 GeV 2 , (26) 

well in agreement with earlier numbers |4S[ . 

The coupling of vector mesons to baryons, formulated 
in the antisymmetric tensor formalism, was discussed in 
Ref. The relation of the generic couplings used in 

that reference to the more conventional vector and tensor 
couplings was given in Ref. |44j . For a more compact 
presentation, we rewrite the terms from Ref. |5l| in a 
SU(2) form, employing immediately standard couplings, 
which results in 



V2 



4m 



My 



{g P [D v , 



N I > 



(27) 



We remark here that we totally neglect p-<j> mixing 
in this analysis (which would be present at least in an 
SU(3) extension). Although the ^-nucleon couplings are 
not as small as the Zweig rule might lead one to expect, 
see e.g. Refs. the mixing is roughly an order of 

magnitude smaller than that between p and cu jU E^l , 
and therefore beyond the accuracy we aim to achieve by 
this model estimate. 



B. Analytical results 

The p-uj mixing contributions to the isospin violating 
Dirac and Pauli form factors can be derived from the 
two diagrams depicted in Fig. 0] They are given in terms 
of the various coupling constants defined in the previous 
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FIG. 4: The two diagrams involving p-uj mixing that con- 
tribute to the isospin breaking form factors. The double lines 
denote vector meson propagators, the cross represents the 
mixing vertex. 



section as follows 
F({t) 
F((t) 



' M v {Mi~ty ' 

' M v {M v -tf ' 

f p m v e ptJ 

(M v - if 
gpKpF^ M v Q pLJ 



(28) 



From these, one can easily derive the leading moments: 



u.d / rp 



dp^p^LU 

(Pi d )rm x = (9^F p ~g p F u )^- 



J pUJ 



(29) 



\P2 ) 



{jJuiKuiFp 9pKp 



2Q P 

Ml 



For the phcnomenological discussion, we will be even 
more interested in the leading moments of the Sachs form 
factors, which are given in terms of the above as 



Pe = Pi 



Am 2 N 



(30) 



u.d u.d , u.d 

Pm = Pi + P2 



In light of the above results, we want to comment on 
the claim made earlier that the saturation of coupling 
constants by vector meson contributions might work bet- 
ter here than for the isospin conserving form factors 
considered in Ref. [44j . The reason is that the addi- 
tional propagator in the mixing case leads to a higher 
power of vector meson masses in the denominators of 
the leading moments, Eq. (|29|l . A heavier pair of isovec- 
tor and isoscalar vector resonances sufficiently close to 
each other in mass to mix, e.g. the p(1450) and the 
w(1420) 1^3, would yield contributions of the same form 
as Eq. 129(1 . but suppressed by a higher power of mass 
ratios My* /My ~ 2 than their unmixed contributions to 
the conventional form factors. 

Finally, we want to comment on the possibility of 
isospin violation other than through mixing. In particu- 
lar, it is possible to construct a mechanism for "direct" 



isospin breaking in the vector-meson-nucleon couplings 
in analogy to Eq. (|12p. 

which results in the p coupling as an isoscalar and the 
ui coupling as an isovector to the nucleons. (Analogous 
terms with the charge instead of the quark mass matrix 
are easily written down.) We disregard this possibility 
for the reason that the vector-meson-nucleon coupling 
strengths are extracted from dispersive analyses on the 
assumption of isospin symmetry (with the exception of 
p-uj mixing in the isovector spectral function), i.e. the 
u> couplings, for instance, are identified as certain pole 
strengths in the isoscalar channel. In this way, an isospin 
breaking p-nucleon coupling would just be taken as part 
of the ui resonance, and vice versa. 



C. Numerical results 

The couplings of the vector mesons to nucleons are 
a rather delicate issue and we prefer to rely as directly 
as possible on data rather than on models. To this 
end, we concentrate on values extracted from disper- 
sive analyses of electromagnetic form factors of the nu- 
cleon, and disregard values extracted from meson ex- 
change models of nuclcon-nuclcon scattering or pion- 
photo-/electroproduction (see e.g. Refs. |5Hl5rll57l l5S jp . 
As it is well known that pure vector meson dominance 
does not yield an adequate description of the isovector 
spectral function, where the two-pion continuum leads 
to a significant enhancement on the left shoulder of the 
p peak |5jj , more recent analyses 0, 0, |(3(], |^| make 
use of the full pion form factor plus irir — * NN partial 
waves. In order to approximately disentangle the spec- 
tral function from Ref. [r3lT | into a non-resonant two-pion 
continuum plus a p contribution, we follow the method 
of Ref. [63] and add a Breit-Wigner parameterization 
of the p resonance to either the chiral one-loop or the 
two-loop |62]| representation of the two-pion cut contri- 
butions. This decomposition is model dependent, but 
probably adequate for a model estimate of low-energy 
constants. The different values in Table |I] give a rather 
consistent picture of the p-nucleon coupling constants on 
an accuracy level of 20-30%. See Appendix f° r the 
relations between various coupling definitions. 

The u> coupling constants are calculated from pure 
zero-width resonance pole residues as found in dispersive 
analyses. The most noteworthy point about the numbers 
in Table [H] is the sign change in in Refs. 0> 113 as 
compared to Ref. [46| and other earlier analyses. While 
the vector coupling g u seems to be determined consis- 
tently (although rather larger than what is inferred from 
AW-scattering |55|,|5(|), the tensor coupling g^K^ is not 
at all, with not even the sign fixed. This uncertainty in 
k w turns out to be by far the dominant uncertainty in 
this analysis. 
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TABLE I: p-nucleon coupling constants and corresponding form factor contributions. "1-loop" and "2-loop" refer to the part 
of the non-resonant two-pion continuum that has been used to extract the p couplings. For details, see main text. 





Ref. [46, 60J 


Ref. [61] + 1-loop 


Ref. [61] + 2-loop 


9p 


4.0 


5.4 


6.2 




6.1 


6.8 


5.1 


10 2 X 


0.9 


1.4 


1.2 


10 2 x -4- [GeV" 2 ] 


0.3 


0.4 


0.3 


10 2 x pi [GeV" 2 ] 


0.3 


0.3 


0.4 


10 2 x 4 [GeV" 2 ] 


3.1 


4.6 


4.0 



TABLE II: tj-nucleon coupling constants and corresponding form factor contributions. For the details, see main text. 





Ref. [46] 


Ref. [47] 


Ref. [63] 


9u, 


41.8 


43.0 


42.2 




-0.16 


0.41 


0.57 


10 2 X K* 


-0.8 


2.2 


3.0 


10 2 x -^r [GeV" 2 ! 


-0.2 


0.6 


0.8 


10 2 x pi [GeV" 2 ] 


8.8 


9.1 


8.9 


10 2 x pi [GeV" 2 ] 


-2.8 


7.5 


10.1 



Considering the leading moments of , we observe 

the following: 

1. K u ' d : The numbers in Tables [jj ITU suggest that, 
as the small tensor coupling of the to is relatively 
enhanced by the larger p decay constant, and 
contribute numbers of similar size. The uncertainty 
in completely dominates the uncertainty of the 
estimate, 

<& = -0-020... + 0.020, (31) 

where we have also included the error in the de- 
termination of the mixing angle in Eq. i|26[l . The 
negative value corresponds to the most recent dis- 
persive analysis jg^. 

2. p^ d : The large vector coupling of the u>, enhanced 
by the p decay constant, completely dominates this 
quantity, such that one has 

(Pe ) mix 3 -{Pe) mix 

Foldy terms hardly play a role, such that even the 
large uncertainty in k u does not spoil the predic- 
tion. We find 

(^) mM = -[ - 07 ---°- 10 ] GeV_2 . ( 32 ) 
where the range includes the uncertainty in Q puJ . 



3. p w ^j : Again, the ui couplings to the nucleon yield 
the largest contributions. The uncertainty is dom- 
inated by the uncertainty in , leading to a range 

(^M)™, = -[0-03...0.15]GeV- 2 . (33) 

The recent analysis |63j| suggests the values larger 
in absolute magnitude within that range. 

D. Discussion 

The above estimates for K u,d , p^ d , p^f correspond to 
low-energy constants that enter the chiral representations 
of the corresponding form factors at 0(p ), 0(p 5 ), 0(p 6 ), 
respectively. As we have only calculated G w ^ uptoO(p 4 ) 
and G^f up to 0(p 5 ), the estimate for n u ' d is the only one 
that is strictly needed for numerical results. However, in 
order to test the potential size of higher-order corrections, 
it may still be very useful to compare the estimates of the 
previous subsection to the chiral loop contributions. 

One downside of the resonance saturation method is 
that, if we want to identify some resonance contribution 
with the finite part of a chiral coupling constant, it is not 
obvious at what scale this identification is to be made. 
It seems to be common wisdom, though, that this scale 
ought to be roughly at the resonance mass; we will there- 
fore choose A = M p . We decompose Eq. (|lb[) according 
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0.04 -| 
0.03 
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-0.01 



LO+NLO, Ktr^ 
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0.1 



0.2 
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CD 
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0.04 -| 
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0.02 -| 
0.01 

o 4 

-0.01 



LO+NLO+VM mixing 
LO+NLO, k ct =0 
LO, K rT =0 







0.1 



0.2 



-t [GeV 



0.3 



FIG. 5: The form factor G u K f(t). The dashed line is the LO 
chiral prediction, the full line includes the NLO order correc- 
tions, both with Kq t (M p ) = 0. The grey band is the chiral 
NLO representation with k°^(M p ) = -0.020 . . . + 0.020. 

to 

K u ' d = K u /{M p ) + K u £{M p ) (34) 

and set Kg T (M p ) — K^ix- Numerically, we find from 
Eq. i|16|) for the chiral part 

K% d {M p ) = 0.025 . (35) 

If we vary the saturation scale, A, in the range 
0.5 ... 1.0 GeV, we find «£> d = 0.020 . . . 0.028, which may 
serve as an indicator of the uncertainty of the resonance 
saturation method as such. These values are all of the 
same order of magnitude as K^ ix in Eq. (|31|l . Even within 
the large error range for K^ d x , however, we predict K u ' d 
to be positive, 

K u ' d = 0.005... 0.045 . (36) 

The most recent values for the coupling constants, with 
a negative K^ ix , lead to a substantial cancellation be- 
tween loop effects and counterterm contributions, and 
altogether to a very small total n u ' d . 

For the magnetic radius term, the two leading chiral 
contributions are unambiguously given in terms of loop 
effects. Evaluating Eq. I|21|l numerically, we find the chi- 
ral prediction at NLO to result in 

(p u ! f) x = (0.05 - 0.04) GeV" 2 = 0.01 GeV" 2 . (37) 

Fig. [S] shows the purely chiral NLO representation of 
G u A f(t) in the range < -t < 0.3 GeV 2 , together 
with the range of counterterm values as estimated from 
Eq. I|31l) . shown as a grey band. For comparison, the LO 
and NLO representations are also shown for Kq^(M p ) = 
0. Although the uncertainty is sizeable, G^f(t) is pre- 
dicted to be positive and smaller than 0.05 in this range. 



FIG. 6: The form factor Gfy (t), notation as in Fig. El The 
grey band now includes the full t-dependence of the mixing 
amplitude. For details, see main text. 

Due to the substantial cancellation between chiral LO 
and NLO contributions to the magnetic radius, the t- 
dependence is very weak. Also, the curvature induced by 
chiral loop effects is minimal. 

Even if the large cancellation between the LO and 
NLO terms for the magnetic radius is considered acci- 
dental, the (formally next-to-next-to-leading) vector me- 
son contribution in Eq. (|33|l is seen to be at least of the 
same order of magnitude, potentially larger than both. 
The vector meson mixing contributions lead to a sign 
change in the radius as compared to the chiral prediction. 
Together with a positive magnetic moment, this means 
the form factors increase in absolute magnitude for non- 
vanishing virtuality in electron scattering experiments, 
t < 0. This is shown in Fig. where, in comparison to 
Fig. [SI we have replaced the low-energy constant contri- 
bution Kq T (M p ) = K^ ix by the full t-dependence of the 
mixing amplitude, Eq. 128|) . As the complete mixing am- 
plitude contains no more parameters than the constant 
at t = 0, the uncertainty band can even get narrower: the 
upper boundary of the band (given by the vector meson 
parameters from Ref. |4^) changes very little and stays 
below 0.05, while the lower boundary (given essentially 
by the parameters from Ref. (6j|) rises with — i. 

We want to emphasize that this combined chiral plus 
vector meson mixing representation as shown in Fig. [B] 
goes beyond strict effective field theory. We believe 
however that the additional i-dependence of the mix- 
ing contributions provides a good estimate of the most 
important higher-order terms that go beyond our chiral 
calculation, for the following reasons: the strongest t- 
dependence from pion loops has to correspond to cuts in 
the low-energy region, and among these, two-pion cuts 
are certainly the most prominent (see e.g. Ref. |59j on 
three-pion cut contributions to nucleon form factors); but 
we have already calculated the two-pion cuts up to NLO, 
and we do not expect higher-order corrections to these to 
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FIG. 7: Contributions to the form factor G u ^ a {t). The full 
line is the LO chiral prediction, the dashed line contains par- 
tial NLO order corrections as discussed in the main text. The 
grey band is the vector meson mixing contribution. 



FIG. 8: The complete form factor G£ d (t). The band com- 
bines uncertainties in the vector meson mixing amplitude as 
well as from higher-order chiral corrections. For details see 
main text. 



be unreasonably large (witness Ref. for two-loop cor- 
rections to the two-pion continuum) . We therefore find it 
reasonable to expect low-energy constants (correspond- 
ing to resonance physics) to dominate the missing pieces 
beyond the NLO chiral representation. We have argued 
earlier that we expect resonance contributions beyond 
p-uj mixing to be minor corrections, and consider their 
possible impact well-covered by the uncertainty bands in 
Figs. El El Note finally that such a combined approach 
was shown to work rather well for the isospin conserving 
form factors in Refs. EH ■ 

We now turn to the electric form factor G E ' (t). It 
can only be predicted unambiguously to leading order 
in ChPT, where one has G u /{t) = F^' d {t). This is 
shown as the dashed line in Fig. [7| A partial higher- 
order correction is given by the chiral contributions to 
t/im 2 N x F^ d {t), see Eq. ©, which is added to the 
leading-order expression for the full line in Fig. [7| and 
can be seen there to be also numerically subleading. 

The leading chiral prediction for the electric radius 
term is, from Eq. (|19l) and Appendix IA II 

<^A = ^^(»'>»). <*> 

Numerically, this amounts to 

( P u /) x = 0.03 GeV" 2 , (39) 

which demonstrates that the vector meson contribution, 
see Eq. (|32|l . is numerically dominant, albeit formally 
subleading. It again leads to a sign change compared to 
the leading chiral prediction. This can also be seen from 
Fig. [7| where the full p—u) mixing amplitude is depicted 
as a grey band, corresponding to the range of vector me- 
son coupling constants yielding the electric radius range 
in Eq. (|32[l . Compared to G^f(t), the band is better 



constrained as the particularly controversial coupling k u 
plays no major role in the electric form factor. Chiral 
and vector meson contributions are combined in Fig. |S1 
where we have taken the partial NLO chiral contribu- 
tions described above as an uncertainty for higher-order 
loop corrections. Errors from both sources were added 
linearly. We consider this band a conservative estimate. 
Due to the bigger mixing amplitude, the total form fac- 
tor is positive, but remains small (< 0.01) and therefore 
well-constrained in the whole momentum transfer range 
considered. 

In order to put these numbers into perspective con- 
cerning the strangeness form factor measurements, we 
want to compare them to some experimental results on 
the latter. Ref. [65| is a recent attempt to combine all 
available world data on parity violating electron scatter- 
ing and perform a best fit for the leading strangeness 
moments. The fit including only leading-order moments 
(e.g. no strange magnetic radius) yields 

k s = 0.12 ±0.55 ±0.07 , 

9 (40) 
p s E = (-0.06 ± 0.41 =f 0.00) GeV , 

while a fit allowing for next-to-leading-order moments 
results in p s M = (0.7 ± 6.8) GeV -2 , i.e. the data are not 
sufficiently accurate yet to pin down the magnetic radius 
to reasonable accuracy. So while the central values of 
Eq. H40|) are already of comparable magnitude to K u ' d , 
p E ,d , the isospin violating moments are still, by a factor 
of 6-10, smaller than the combined uncertainties on the 
strangeness moments. 

In Table ITTT1 we also compare to a few selected indi- 
vidual experimental numbers on strangeness form fac- 
tors P, 0, 01 . We contrast those results with bands for 
the isospin violating form factors in the same kinematics, 
i.e. for the same combination of electric and magnetic, 
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TABLE III: Comparison of selected experimental measurements of strange form factors to the results of this work. The column 
"electric/magnetic" shows the linear combination of strange electric and magnetic form factors as well as the momentum 
transfer in the respective kinematics, the column "G s " the corresponding result. The column u Q u > d " gives the range of the 
isospin breaking form factor contributions for the same linear combination, in the same kinematics, as predicted in this work. 

experiment electric/magnetic G s G u ' d 

SAMPLER G A / (-0.1 GeV 2 ) 0.37 ± 0.20 ± 0.26 ± 0.07 0.02 ... 0.05 

A4 [3] [G E +0.106 G M ] (-0.108 GeV 2 ) 0.071 ± 0.036 0.004 ... 0.010 

HAPPEX [5] [G E +0.080 G M ] (-0.099 GeV 2 ) 0.030 + 0.025 + 0.006 + 0.012 0.004 ... 0.009 



and the same momentum transfer t. The error band cor- 
responds to a worst-case combination of the bands dis- 
cussed in connection with Figs. H3 El The conclusion is 
similar here: the uncertainty from isospin breaking is still 
smaller than the overall experimental error, but is rele- 
vant as measurements become more and more precise. 

V. SUMMARY AND CONCLUSIONS 

In this paper, we have re-investigated the isospin 
breaking vector form factors of the nucleon that are in- 
creasingly becoming a necessary ingredient to the precise 
extractions of the strange vector form factors. Our find- 
ings can be summarized as follows: 

1. We calculate the isospin breaking electric and mag- 
netic form factors G u ^j M (t) to leading and next-to- 
leading order in two-flavor chiral perturbation the- 
ory, correcting some errors in Ref. |29|. Up to this 
order of accuracy, all loop effects are due to the 
proton-neutron mass difference. Only one combi- 
nation of unknown low-energy constants enters the 
chiral representation of these form factors, namely 
a contribution to the isospin breaking magnetic mo- 
ment. Leading and next-to-leading-order contribu- 
tions to the magnetic radius are seen to cancel to 
a large extent, leading to a very weak momentum 
dependence. 

2. In order to estimate the missing counterterms, we 
employ the method of resonance saturation. We 
provide analytic expressions for the isospin break- 
ing magnetic moment as well as for (formally sub- 
leading) electric and magnetic radii in terms of var- 
ious vector meson coupling constants and the p—u) 
mixing angle. While most of these phenomenolog- 
ical parameters are known to apt precision, by far 
the most uncertain coupling is the w-nucleon ten- 
sor coupling, which dominates the rather broad er- 
ror band in our prediction for the isospin breaking 
magnetic form factor. 

3. For the electric and magnetic radii where counter- 
terms only come in at subleading orders, we find 
that the combination of the relatively small nucleon 



mass difference, which is responsible for the chiral 
loop contributions, and the strong vector-meson- 
nucleon couplings, which numerically enhance the 
effect of low-energy constants, tends to upset the 
hierarchy suggested by chiral power counting. The 
vector meson mixing contributions lead to sign 
changes in the radii compared to the purely chiral 
predictions. 

4. We combine chiral loop contributions and the full 
vector meson mixing amplitudes in a phenomeno- 
logical approach, adding the various uncertainties 
to produce conservative error bands. We find that 
^E/M are both positive in the momentum trans- 
fer range < — t < 0.3 GeV 2 , with approximate 
upper limits G u A f{t) < 0.05, G u /{t) < 0.01. 

In order to sharpen our findings and make the prediction 
for the isospin breaking vector form factors even more 
stringent, the predominant task would be to improve on 
our knowledge of the w-nucleon tensor coupling, if feasi- 
ble. Higher-order chiral calculations would be extremely 
tedious as they would involve two-loop diagrams with pi- 
ons and photons, and are unlikely to be of high numerical 
relevance compared to the resonance physics considered 
in this article. 
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APPENDIX A: CHIRAL BACKGROUND 
1. Loop functions 

Here we spell out the explicit forms of the loop func- 
tions used in Section UlI CI 



7o (t) 
73 (t) 



i 



2M. 



arctan 



16tt V V-t 
I / 4Ai 2 _ t 



16?r V 4M„ 
1 



2M, 
arctan 



16tH 



(7 log 



t 

g + 1 

(T- 1 



2M„ 
2 



(Al) 



where cr = yT— 4M%/t. 7o/aM are linked to the loop 
functions 70/3 (*) used in Ref. |36| via 

7b(«) = ^(7o(t)-7o(0)) , 7s(i) = M~ 1 (73 (*) — 73 (0)) . 
The leading ^-dependence of these functions is given by 



7o (*) 
73 W 



192ttM2 
t 

9671-2M 2 4 



0(i 2 ) , 
+ 0(i 2 ), 
0(i 2 ) . 



(A2) 



We note that the functions in Eq. (|A1|) are the heavy- 
baryon loop functions that were reproduced from the in- 
frared regularized loops by strict expansion in chiral pow- 
ers. As we are only interested in the form factors in the 
space-like region, we disregard the complications ensuing 
from the anomalous threshold of the triangle diagram. 
We want to emphasize, though, that these representa- 
tions do not reproduce the correct threshold behavior of 
the spectral functions. 



2. Separate diagrams 



For completeness, we show here the various contribu- 
tions to Ff, 2 (Q), if , 2 (0) from the diagrams displayed in 

Fig. [J Of course if (0) = if (0) = in the sum of all 
diagrams is a necessary requirement. The pref actor C is 
as defined in Eq. (fHfl . 



if (0) [(a)] =-2if [(b)] =-2if [(c) 

#)]=-#)]~c, 



m N 



C , 



(A3) 



if (0) [(a)] 

^ [(d)] 
if [(d)] 



l + 21og- 



m N 



C , 



m N 

~3 F ^^= ^ C ' 



1 



J* [(e)] 



w(k p + K n )M w 



m N 



C . (A4) 



APPENDIX B: 



VECTOR MESON COUPLING 
CONSTANTS 



In this appendix, we briefly spell out how to relate 
the vector-meson-nucleon coup ling constants used in this 
paper to those in Refs. I471 1621 . 

In Ref. [§2 , the strength of the p contributions to the 
isovector spectral functions of the Sachs form factors is 
parameterized in terms of two coupling constants &e,m 
that are given in terms of g p , n p according to 



9 P Fp 
2M 



4m 2 v 



9 P Fp 
2M 



(1 



(Bl) 



In Ref. 62], the numerical values bs = 1.0, 6m = 3.6 were 
extracted. Using th e up dated empirical isovector spectral 
function from Ref. [61| , we find the slightly shifted num- 
bers = 1.1, 6m = 3.7, while fitting the phenomeno- 
logical p contribution with just the one-loop two-pion 
continuum leads to 6g = 1.2, 6m = 4.2. These numbers, 
together with Eq. {BTJ, feed into Table |T] 

The lu coupling constants are given as pole residues 
2 m the isoscalar spectral functions of Dirac and Pauli 



form factors in Refs. |46l 1471. 63] . The coupling constants 
g u , Kaj can be calculated from these as 



2a% 

f~mL 



a? 



(B2) 



The pole residues from Ref. [6J| used for Table HU are 
af = 0.752 GeV 2 , a% = 0.425 GeV 2 . 
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